Ground states with cluster structures in a frustrated Heisenberg chain 
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We examine the ground state of a Heisenberg model with arbitrary spin S on a one-dimensional lattice 
composed of diamond-shaped units. A unit includes two types of antiferromagnetic exchange interactions 
which frustrate each other. The system undergoes phase changes when the ratio A between the exchange 
parameters varies. In some phases, strong frustration leads to larger local structures or clusters of spins 
than a dimer. We prove for arbitrary S that there exists a phase with four-spin cluster states, which was 
previously found numerically for a special value of A in the S — 1/2 case. For S — 1/2 we show that there 
are three ground state phases and determine their boundaries. 
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PACS numbers: 05.30.-d, 75.10.Jm, 75.50.Gg 

Effects of frustration in quantum antiferromagnets are 
of great current interest in solid state physics. In a clas- 
sical system, strong frustration obstructs a simple anti- 
ferromagnetic (AF) ordering and produces another mag- 
netic order characteristic of each system, e.g. the 120° 
structure or a spiral Jl[] . In a quantum system, the inter- 
play of quantum fluctuations and frustration makes the 
situation more complicated. There may appear an exotic 
ground state which does not have magnetic order and has 
no classical analog. A typical example is the complete 
dimer state in the Majumdar- Ghosh (MG) model g|. 

There are two types of frustrated quantum spin sys- 
tems. One is a system the classical version of which has 
a locally stable spin configuration in the ground state; 
i.e. any local deformation for a spin configuration always 
raises the energy. We say that such a system is elastic. 
The MG model, the model with linearly decreasing AF 
interactions HQ and the AF Heisenberg model on the 
triangular lattice are of this type. For the other type, the 
classical version of a system has ground-state spin con- 
figurations which can be locally deformed without rais- 
ing the energy. Then the set of these configurations is 
a manifold with dimensions proportional to the system 
size. We say that a system of this type is floppy. The AF 
Heisenberg models on the A chain g, the double chain 
with diagonal interaction jfj) and the kagome lattice |t]] 
are of this type. Some other floppy spin systems are also 
seen in ||. 

Specially interesting is floppy systems. For example, 
it is argued that the kagome antiferromagnet has a mys- 
terious peak in a low temperature part of the specific 
heat 0. In spite of quite a few theoretical studies on 
this system, its quantum ground state and the low tem- 
perature thermodynamic properties are hardly clarified. 
The difficulty of the problem originates from the floppi- 
ness. In a classical floppy system, the configuration of a 



local set of spins can be deformed without affecting the 
other part. In the corresponding quantum system the 
local set may form a nearly closed state, or a cluster. A 
part surrounded by such clusters forms another cluster. 
Thus the total wave function becomes approximately of 
a direct product form. To confirm this picture, it is im- 
portant to examine a simple floppy model in which the 
structure of the ground state is clearly seen. 

The singlet dimer is the smallest cluster. Ground 
states with dimer structures are possible in both elastic 
and floppy systems with S = 1/2. The difference appears 
in their low-lying energy spectra. An elastic quantum 
system has usually unique or finitely degenerate ground 
states. The lowest excitation mode is expected to have an 
energy gap and a dispersion both of the order of the typ- 
ical exchange parameters. On the other hand, unusual 
situations such as a macroscopic degeneracy of the dimer 
ground state || or dispersionless low energy excitations 
H occur in floppy systems. These seem to be two differ- 
ent quantum manifestations of the classical floppiness. 

In this paper we report exact ground states with larger 
clusters than dimers in a frustrated AF chain. The sys- 
tem contains a parameter which regulates the strength of 
the frustration, and can be either elastic or floppy accord- 
ing to the value of the parameter. Wave functions with 
periodic cluster structures appear in the floppy regime. 

The system is described by the Hamiltonian: 
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hi = J(Si + S l+1 ) ■ (Tf } + TH + J'T] 1 ' ■ T\*> (lb) 

where Si, Tj and are spins with magnitude S in 
the i-th unit cell. We assume that J is positive. The 
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sub-Hamiltonian hi corresponds to a diamond-shape and 
so we call the total lattice shown in figure 1(a) simply a 
diamond chain. 

The system was studied earlier briefly by Sutherland 
and Shastry as an example of the superstability || . The 
case of 5=1/2 and J' = 2 J was discussed as an example 
in the context of a general method to derive spin models 
with complete dimer ground states The cases with 
J' = ± J were studied by Long et al 0. They numerically 
found the ground state with four-site clusters for 5 = 1/2 
and J' = J. 

Recently it was reported that the magnetic properties 
of a mixed molecular crystal of organic mono- and di- 
radicals are described by the Hamiltonian (l)for5 = 1/2 
with J ~ 30K and J' ~ -20K @. 

We treat the whole region of J', though our main in- 
terest is in the frustrated region where J' is positive. For 
a given 5 the ratio A = J'/ J is the only parameter which 
determines the properties of the system. We assume pe- 
riodic boundary conditions in the following, though most 
of the results do not depend on this assumption. 

Before proceeding to the full quantum treatment, we 
examine the ground state spin configuration in the clas- 
sical limit (5 — ► oo with J5 2 finite). The ground state is 
given by a spin configuration which minimizes the energy 
of all the unit diamonds. Equation (lb) is rewritten as 

ht = + S,) 2 + (A - 1)T 2 - S 2 - 2A5(5 + 1)} (2) 

by using T, = t| 1} + T ( (2) and Si = Si + S l+1 . In the 
classical limit 5(5+1) is replaced by 5 2 in (g). For A < 1 
it is seen that the energy is minimized when the following 
conditions are satisfied: |Tj| = 25, |Sj| = 25 and |Tj + 
Sj| =0. These conditions imply that T^ 1 - 1 and T[ , and 
Sj and Si + i are parallel to each other, respectively, and 
further (j=l or 2) and Sj are antiparallel to each 
other. So this is a ferrimagnetic state where all the Si's 
(Tp^'s) are aligned parallel (antiparallel) to an axis. The 
ground state energy and the magnetization are given by 
— (4 — A) JS 2 and 5 per unit diamond, respectively. The 
ground state is elastic in this region. 

For A > 1, another expression of hi is useful, i.e. 
hi = (JA/2){(T. 1 + S l /A) 2 -(S l /A) 2 -25(5+l)}. The en- 
ergy minimum is realized if |S<| = 25 and |Tj+Sj/A| = 0. 
The former condition leads to all the Si's being parallel 
to an axis. The latter condition is satisfied if and only 
if A > 1, and means that and T^ 2 ' form a trian- 
gle with Si/A. They make an angle 8 — cos^ 1 1/A to the 
axis through Si and may be rotated about this axis simul- 
taneously without raising the energy. The ground state 
contains N degrees of freedom of free rotations and hence 
is floppy. The ground state energy is — (A + 2/ A) JS 2 per 
unit diamond. The ground state is again ferrimagnetic 
except for A = 2 with a magnetization |2/A— 1|5 per unit. 



Thus in the classical limit the diamond chain has two dif- 
ferent ground state phases, the elastic one for A < 1 and 
the floppy one for A > 1. 

We turn to the quantum case with general 5. Hereafter 
we take an energy unit of J = 1 . We find easily from (|^) 
that [T 2 , H] = 0; i.e., for all i, Tf = Ti{Ti + 1) are good 
quantum numbers (Ti = 0, 1, ... , 25). With fixed {Ti}, 
the original problem of 37V spins reduces to a problem of 
a linear chain with 27V spins, where the (2i — l)-th site is 
occupied by the spin Si and the 2i-th by Tj . The energy 
of J'-bonds (i.e. a part of the energy proportional to A) 
is determined solely by {Ti}. It should be noted that if 
Ti = then there is no interaction between the left and 
the right side of Ti and the whole lattice is decoupled. 

Let us first consider the eigenstates of an isolated unit 
diamond described by hi. In the lowest eigenstate of hi 
for a given Ti, the energy is 

Ei(Ti) = Til-fa + 1) - (25 + 1)] - A5(5 + 1) (3) 

and the total spin is 25 — T . Then we obtain, 

Lemma 1. For T > 1 and for A > Ad(T) where 
A D (T) = 2(25 + 1)/(T + 1), it holds that T L ^T for any 
tin the ground state of the diamond chain. 

proof The total Hamiltonian is divided as H = H'+ 
h m , where H' is the sum of h^s with i ^ m. Let us 
take a state |\&o) whose wave function is given by the 
direct product of the ground state wave function of H' 
and the singlet wave function of T m . Then (^ol^l^o) = 
Ei(0) + E', where E' is the ground state energy of H' . 
Let |*) be any state with T m = T > 1. Then (^\H\^) = 
(9\hm\9) + (*\H'\V). Clearly (*|/i m |#) > Ey{T) and 
(tfj-ff'l*) > E'\ also Ei(T) > £i(0) for A > A D (T). 
Therefore {9\H\^B) > (^q\H\^ ) and \^>) can not be the 
ground state of H . QED 

Since Ad(T) decreases with T from Ad(1) = 25 + 1 to 
Ad (25) = 2, we obtain the following result. 

Proposition 1. For A > 25 + 1, T t = for alH in 
the ground state. That is, all pairs of and form 
singlet dimers. 

In this state all Si's are decoupled from other spins 
and behave as free spins. So we call this state a dimer- 
monomer (DM) state (figure 1(b)). Due to the free spins 
there is a (25 + 1) -fold degeneracy in the DM state. 

For A < 25+1, the DM state is not a ground state 
since one can lower the energy by introducing an isolated 
T = 1 in the DM state. For 5 = 1/2 the ground state is 
composed of Tj's with their magnitude or 1. It is also 
true for 5 > 1 if 25 + 1 > A > A D (2) = 2(25 + l)/3. 
For example, let us assume the configuration of {Ti} 
as {1010111011110110} in the chain with 16 diamonds. 
Since zero Ti's decouple the system to n-diamond clus- 
ters, the lowest energy for this configuration is simply 
given by 5£i(0) + 2E X + E 2 + E 3 + E 4 , where E n (n>\) 
is the lowest energy of a cluster with n diamonds for the 
configuration with all Ti's equal to 1. Taking account of 
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an i?i(0) accompanying each n-diamond cluster, the av- 
erage energy of this cluster per diamond measured from 
the DM state is given by e„ = (E n - nEi(0))/(n + 1). 
If the minimum value of e„ occurs at n = n m , then the 
ground state in the thermodynamic limit (N — * oo) is 
realized by a regular array of n m -clusters with isolated 
T = O's between them. A finite N may cause a mis- 
match of the cluster size. If the minimum e n is realized 
at more than one value of n, then the ground state will 
have a macroscopic degeneracy. It seems plausible that 
the ground state is the state with 1-diamond clusters on 
every other diamonds for A less than but close to 25 + 1, 
since a simple variational argument leads to E n > nE\ 
and usually the difference between both sides of the in- 
equality is fairly large. The wave function of this state 
is a direct product of those of four spin clusters (T = 1) 
and dimers (T = 0) as shown in figure 1(c) || . Therefore 
we call this state a tetramer-dimer (TD) state fTl|| . 

Proposition 2 There exists a finite region Aci < 
A < 25 + 1 in which the TD state is the ground state of 
the diamond chain in the thermodynamic limit. 

proof We show that e n is the minimum only at 
n=l for a region Aci < A < 2S + 1. We denote the 
ground state energy of a nearest neighbor AF linear chain 
with n+1 spins of magnitude 5 on odd sites and n spins 
of magnitude unity on even sites as E n . Then we have 
E n — nE\ (0) = E n + nX. Let us first assume 5 = E 2 — 
2E\ > 0. If n is even, the Hamiltonian for E n is divided 
into n/2 sub-Hamiltonians, each equivalent to that for 
E2 ■ A variational argument gives an upper bound on E n 
as nE\ + n5/2 for even n, and as nE\ + (n — l)S/2 for 
odd n (> 3). Since e\ = [A — (2S + l)]/2, a lower bound 
on e„ — ei is given as 2 (n+i) — (^5 + -0 + wmcri ls 
positive for A > 25 + 1 — 6 and n > 2. So we obtain 
Aci < 25 + 1 - S. 

Next we assume that 6 = and deduce a contradiction. 
We divide the Hamiltonian for E2 with five spins as hi + 
h%, where each one is the Hamiltonian for three spins, 
two of magnitude 5 and one of magnitude unity in the 
middle. A simple variational argument proves that d > 
and the assumption implies that the ground state of 
hi + hi is also the ground state of hi. The Lieb-Mattis 
theorem holds in this case and implies that the total 
spin of the ground state is |35 — 2\ and that for hi is 
25 — 1. The theorem also tells us that the ground state 
with the z-component of the total spin 5^ contains all 
the 5 z -diagonal states compatible with S^. The state 
5, 0, 5, —1, 5 — 1) is contained in the ground state with 
$f = 35 — 2. But the fact that the total spin of three 
spins on the left must be 25 — 1 eliminates the presence of 
the above state in the ground state and proves the failure 
of the assumption. QED 

Remark We can show that Aci < 1 for 5 = 1/2. 
The proof will be given elsewhere |l3| . 

For each 5, we need numerical studies of finite size 



systems to determine the precise value of Aci and to find 
the ground state structure for A < Aci- The analysis for 
5 = 1/2 and 1 will be given later. 

For negative or positive but small A, we expect that the 
ground state is a ferrimagnetic state corresponding to the 
complete ferrimagnetic ground state in the classical limit. 
For A < we have a general result. 

Proposition 3 For A < the ground state is fer- 
rimagnetic, i.e. it has both ferromagnetic and AF long- 
range order. For all i, T = 25 in the ground state. 

proof We divide the total lattice into two sublat- 

( 7*1 

tices, where Si's are on the A sublattice and T s on 
the B sublattice. For A < the Lieb-Mattis theorem |13] 
implies that the total spin of the ground state is given 
by NS; i. e. ferromagnetic long range order (FLRO) ex- 
ists. The positive (or negative) definiteness of the ground 
state wave function in terms of the 5 2 -diagonal basis im- 
plies that AFLRO is not less than FLRO Q and that 
T = 25 for all i. QED 

The ferrimagnetic ground state phase is expected to 
extend to positive A and we determine the phase bound- 
ary for 5 = 1/2 by numerical methods. 

There is macroscopic degeneracy of the ground states 
at A = 25+1. Since Ei(l) = Ei(0) at this value 
of A, any configuration of {T} with only T = 1 or 
and satisfying {Ti,T + i} ^ {1, 1} for any i gives ground 
states. We estimate the total degeneracy D by using an 
asymptotic expansion for large N as log D ps N log[(25 + 
l)'iyi + 4/a + l)/2] where a = (25 + 1)7(45 - 1). 
The ratio of the residual to the total entropy is 0.424 
for 5 = 1/2 and decreases with 5 to 1/3 for 5 = 00. 
Macroscopic degeneracy of the ground state can occur at 
other phase boundaries. 

Let us inspect the case of 5 = 1/2. We estimate the 
minimum energy of a cluster composed of n unit dia- 
monds where all TVs are unity. For this purpose, we 
numerically calculated the ground state energy E n of fi- 
nite linear chains with n + 1 spin one-halves and n spin 
unities alternatingly aligned. By employing the Lanczos 
technique we obtained the energy up to n = 7. The ob- 
tained data of E n fit nicely to an assumed asymptotic 
formula E n / (n + 1) ~ + a/(n + 1) with = —1.454. 
We have estimated the values of A where the state with n- 
diamond clusters gives the same energy with the TD state 
by using the obtained E n 's. The values are 0.763012, 
0.819171, 0.847024, 0.862127, 0.871321, and 0.877464 for 
11 = 2, 3, . . . , 7. From the monotonically increasing be- 
havior of these values, and the linear dependence of e n 
on A we conclude that the TD phase changes directly to 
the ferrimagnetic phase where T = 1 for all i. The crit- 
ical A is estimated as 0.909 from the value of iao. Thus 
we have completely determined the ground state phase 
diagram of the diamond chain for 5 = 1/2. It consists of 
three phases; the DM phase for A > 2, the TD phase for 
0.909 < A < 2 and the ferrimagnetic phase for A < 0.909 
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(figure 2). 

In the case of S = 1 we have obtained ground states 
with clusters larger than tetramers. The ground state 
changes its character successively from the DM state to 
the TD state at A = 3, then to the heptamer (n=2)-dimer 
state at A = 2.660, to the state with 3-diamond clusters 
at A = 2.583 and suddenly to the state where 7$ = 1 for 
all i (i.e. n = oo) at A = 2.577. Therefore the ground 
state phase diagram has at least six phases. The detail 
will be reported elsewhere Jl3| . 

Above we have shown that the spin system on the di- 
amond chain have several ground state phases with clus- 
ter structures There appear clusters larger than 
conventional dimers. For S — 1/2 the TD state is es- 
pecially interesting since it is nonmagnetic. We believe 
that these cluster structures are quantum manifestations 
of the classical floppiness. The results might give some 
insight into the ground state of other floppy systems such 
as the kagome antiferromagnet. 
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Figure Captions 

Figure 1. (a) The diamond chain. The circle represents 
a spin with magnitude S and the solid (dashed) line 
represents exchange parameter J(J'). (b) The dimer- 
monomer (DM) state. The unshaded oval represents a 
dimer (Tj = 0). There are free spins (monomers) on 
the sites not closed by ovals, (c) The tetramer-dimer 
(TM) state. The shaded oval represents a triplet pair 
(Ti = 1) and the closed loop including four spins repre- 
sents a tetramer. 

Figure 2. Phase diagram for S — 1/2 in the parameter 
space of A (= J' /J). 
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